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Abstract 

Analytical and geometrical properties of generalized power-law 
(GPL) density profiles are reviewed, and special effort is devoted to 
the special cases where GPL density profiles reduce to (i) a double 
power-law (DPL), and (ii) a single power-law (SPL). Then GPL den- 
sity profiles are compared with simulated dark haloes (SDH) density 
profiles, and nonlinear least-squares fits are prescribed, involving five 
parameters (a scaling radius, tq, a scaling density, po, and three ex- 
ponents, a, P, 7), which specify the fitting profile (RFSM5 method). 
More specifically, the validity of a necessary condition for the occur- 
rence of an extremal point, is related to the existence of an intersection 
between three surfaces in a 3-dimension space. Using the algorithm 
makes also establish that the extremal point is a fiducial minimum, 
while the explicit calculation of the Hessian determinant is avoided to 
gain in simplicity. In absence of a rigorous proof, the fiducial minimum 
can be considered as nothing but a fiducial absolute minimum. An 
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application is made to a sample of 17 SDHs on the scale of cluster of 
galaxies, within a flat ACDM cosmological model (Rasia et al. 2004). 
In dealing with the averaged SDH density profile (ADP), a virial ra- 
dius, r^r, equal to the mean over the whole sample, is assigned, which 
allows the calculation of the remaining parameters. The following re- 
sults are found, (i) A necessary condition for the occurrence of an 
extremal point is satisfied for 8 sample haloes, and is not for the re- 
maining 9 together with ADP. In the former alternative, an extremal 
minimum point (EMP) may safely exist. In the latter alternative, 
the occurrence of an EMP cannot be excluded, but only a non ex- 
tremal minimum can be determined, (ii) The occurrence of an EMP 
implies a sum of square residuals which is systematically lower than 
its counterpart deduced by use of numerical RFSM5 methods. Ac- 
cordingly, EMPs may safely be thought of as absolute minima. With 
regard to sample haloes where no EMP is detected, the above result 
maintains in 4 cases (including ADP), while the contrary holds for 
the remaining 5 cases, (iii) The best fit (with no EMP detected) is 
provided by DPL density profiles for 3 sample haloes. In addition, 
DPL density profiles make a rough, but viable approximation in fit- 
ting SDH density profiles. The contrary holds for SPL density pro- 
files. No evident correlation is found between SDH dynamical state 
(relaxed or merging) and asymptotic inner slope of the logarithmic 
density profile or (for SDH comparable virial masses) scaled radius. 
Mean values and standard deviations of some parameters are calcu- 
lated and, in particular, the decimal logarithm of the scaled radius, 
^yir, reads < log^^^ >~ 0.59 and 0.59, the standard de- 

viation exceeding by a factor 3.3-3.4 its counterpart evaluated in an 
earlier attempt using NFW (Navarro et al. 1995, 1996, 1997) density 
profiles (Bullock et al. 2001). If a large dispersion still maintains 
for richer samples, in dealing with analytical RSFM5 methods, a low 
dispersion found in A?^-body simulations seems to be an artefact, due 
to the assumption of NFW (or any equivalent choice) density profile. 
It provides additional support to the idea, that NFW density profiles 
may be considered as a convenient way to parametrize SDH density 
profiles, without implying that it necessarily produces the best possi- 
ble fit (Bullock et al. 2001). With regard to RFSM5 methods formu- 
lated in the current paper, the exponents of both the best fitting GPL 
density profile to ADP, (a,/?, 7) « (0.3,4.5, 1.5), and related averages 
calculated over the whole halo sample, {a,P,j) (3.0,3.9,1.3), are 
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far from their NFW counterparts, (a,/3,7) = (1,3, 1). The last result, 
together with the large value of the standard deviation, (Jsiog^yiri is 
interpreted as due to a certain degree of degeneracy in fitting GPL to 
SDH density profiles. If it is a real feature of the problem, or it could 
be reduced by the next generation of high-resolution simulations, still 
remains an open question. Values of asymptotic inner slope of fitting 
logarithmic density profiles, are consistent with results from recent 
high-resolution simulations (Dicmand ct al. 2004; Rccd ct al. 2004). 
keywords - cosmology: dark matter - galaxies: clusters: general. 

1 Introduction 

Recent observations from anisotropies in the cosmic microwave background, 
large-scale structure surveys, Hubble parameter determinations, and Type la 
supernovae results, allow narrow ranges for the values of cosmological quanti- 
ties (e.g., Sievers et al. 2003; Rubino-Martin et al. 2003; Spergel et al. 2003). 
The related (ACDM) cosmological model is consistent with a bottom-up pic- 
ture (hierarchical clustering) of dark matter haloes, where smaller systems 
formed first from initial density perturbations and then merged with each 
other to become larger systems, or were tidally disrupted and accreted from 
previously formed larger systems. 

A two-parameter functional form for the halo profile, where a scaled den- 
sity depends on a scaled radius and three exponents are specified, has been 
proposed by several authors (e.g., Navarro et al. 1995, 1996, 1997, hereafter 
quoted together as NFW; Moore ct al. 1998, 1999; Rasia et al. 2004). Gen- 
erally speaking, it may be considered as a special case of the 5-parameter 
family (Hernquist 1990): 



where po and tq are a scaling density and a scaling radius, respectively, and 
all the exponents are positive in the case under discussion. 

The density profile, expressed by Eq. (1), reduces to a power-law both 
towards the centre, r — » 0, and towards infinite, r — >■ -|-oo, where the expo- 
nent equals —7 and — /5, respectively. It may be conceived as a generalized 
power-law and, in the following, it shall be quoted as GPL density profile. 
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On the other hand, matter distribution within a simulated dark matter halo, 
in the following, shall be quoted as SDH density profile. 

For fixed exponents, one among the two remaining free parameters, the 
scaling density and the scaling radius, may be related to the mass and the 
radius of the virialized region. For further details see e.g., NFW; Bullock et 
al. (2001); Rasia et al. (2004); Caimmi et al. (2004; hereafter quoted as 
CMV). 

Strictly speaking, GPL density profiles cannot be considered as universal, 
in the sense that a fixed choice of exponents, (a,/?, 7), does not provide 
the best fit to an assigned set of SDH density profiles (e.g.. Bullock et al. 
2001; CMV). A weaker conception of universal, GPL density profiles, might 
be related to the insensitivity of the shape to halo masses, initial density 
fluctuation spectra, or cosmological parameters (NFW; Hess et al. 1999; 
Klypin et al. 2001; Bullock et al. 2001; Fukushige & Makino 2001, 2003; 
Zhao et al. 2003). 

However, no general consensus still exists on this point. In particular, 
some authors found an inner slope which depends on the power spectrum of 
the initial density perturbation (Syer & White 1998), or on the mass (Jing & 
Suto 2000; Hiotelis 2003; Ricotti 2003). Further support to the above results 
is provided by recent, high-resolution, A^-body simulations (Fukushige et al. 
2004; Navarro et al. 2004). In addition, a certain degree of degeneracy 
occurs, with regard to the exponents, (a,/?, 7), in fitting various GPL to 
SDH density profiles, in the whole range of resolved scales (e.g., Klypin et 
al. 2001). For a more detailed discussion, see CMV. 

According to the above remarks, GPL density profiles may be considered 
as a convenient way to parametrize SDH density profiles, without imply- 
ing that it necessarily produces the best possible fit (Bullock et al. 2001). 
The search of the best fitting GPL to (one or more) SDH density profiles, 
makes the aim of the current paper. To this respect, different fitting cri- 
teria exist, such as minimizing the maximum fractional deviations of the 
fit, max|log(pGPL/Pfe) - ^og{psDH/ Ph)\ (e.g., Klypin et al. 2001); the sum 
of the squares of absolute^logarithmic residuals, = Z][log(pGPL/p/i) — 
\og{psDH / Ph)Y (e-g-, Bullock et al. 2001); the sum of squares of relative 

^The term "absolute" here has not to be intended as "absolute value" , but as opposite 
to "relative". More precisely, yi — y{xi) is an absolute residual, while \yi — y{xi)]/y{xi) is 
the corresponding relative residual. 
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residuals, Y\{psDH/Ph - Pgpl/ Ph) / {pgpl/ Ph)f (e.g., Fukushige & Makino 
2003, 2004), where ph is a normalization value. For a more detailed discus- 
sion see e.g., Tasitsiomi et al. (2004). 

The hierarchical collapse of dark matter into virialized haloes is likely to 
have played a key role in the formation of large-scale objects, such as galaxies 
and clusters of galaxies. The halo density distribution has a direct dynamical 
role in determining the observable parameters of the baryonic subsystems. 
Then further investigation even on the best fitting GPL to SDH density 
profile, appears to be important. 

To this aim, a nonlinear least-squares method is used in the current paper. 
The related, general procedure, shall be quoted in the following as RFSM5 
(Residuals Functions Sum Minimization within a 5-dimension hyperspace) 
method. The radius of the virialized region is taken from computer outputs, 
for selected cosmological parameters, which allows to deduce the mass of the 
virialized region, and the remaining quantities of interest. Further details 
may be found in an earlier attempt (CMV), where fiducial minima of the 
sum of absolute values and squares of residuals are calculated for assigned 
SDH density profiles, using a 5-dimension hypergrid within a 5-dimension 
hyperspace. 

An analytical procedure, instead of the above mentioned numerical, shall 
be exploited in the current paper. The problem may be formulated from both 
an analytical and geometrical point of view. The sum of square residuals is 
a function of 5 variables in a 6-dimension hyperspace, w — F{ro, po, a, (3, 7) 
which, in turn, represents a 5-dimension hypersurface. The extremal point of 
absolute minimum defines the best fitting GPL to SDH density profile under 
consideration. If more than one extremal points of absolute minimum exist, 
then degeneracy occurs. For further details, see CMV. 

According to a theorem of mathematical analysis, a necessary condition 
for the occurrence of an extremal point of minimum, is a null (or undefined) 
value of all the first partial derivatives of an assigned function, evaluated 
at that point. A study of the Hessian determinant enables to establish ad- 
ditional features (minimum, maximum, hyperflexion, hypersaddle) . Finally, 
one has to ensure that the minimum is an absolute minimum. For further 
details on the analytic procedure see e.g., Smirnov (1969, Ghap. V, §§ 2.2-2.6). 

In fitting GPL to SDH density profiles, our attention shall be restricted 
to the sum of square residuals, as absolute values exhibit a discontinuity of 
the first derivative for null arguments, which makes a solution of the problem 
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much more difficult. An algorithm will be formulated and used, for searching 
points where the above mentioned necessary condition for the occurrence of 
an extremal minimum, is satisfied. Then it will be inferred, by comparison 
with numerical results found in an earlier attempt (CMV), that the related 
point is a fiducial minimum. On the other hand, no proof will be produced, 
that the point under discussion is an absolute minimum. This, and other 
improvements, are left to mathematicians. 

The fit of GPL to SDH density profiles, which has deeply been inves- 
tigated in CMV, is reviewed in section 2. The analytical RFSM5 method 
used for fulfilling the necessary condition for the occurrence of an extremal 
point of minimum, is presented in section 3 and derived in Appendix A. The 
subject of section 4 is an application to a sample of 17 SDHs and the related 
mean density profile, on the scale of clusters of galaxies, taken from Rasia et 
al. (2004), and already used in CMV. Then the results are discussed. Some 
concluding remarks are drawn in section 5. 



2 Fitting GPL to SDH density profiles 

A detailed investigation on GPL and SDH density profiles, together with their 
comparison, has been performed in an earlier attempt (CMV), to which an 
interested reader is addressed. Only what is relevant for the current paper, 
or has not been mentioned in CMV, shall be reported here. 



2.1 GPL density profiles 

Plotting GPL density profiles on a logarithmic plane, (0 log^ log/), neces- 
sarily implies use of dimensionless coordinates, defined as: 

e = /(0 = f; (2) 

where ^ and / can be conceived as a scaled radius and a scaled density, 
respectively. Accordingly, Eq. (1) reduces to: 

independent of the scaling parameters. 
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As Eqs. (2) and (3) imply null density at infinite radius, the mass distri- 
bution has necessarily to be ended at an assigned isopycnic surface, which 
defines a truncation radius, R. The mass within the truncation isopycnic 
surface is (Caimmi & Marmo 2003; CMV): 

4:71 

M = M{R) = -jr-oPol^mas = Mommas ] (4) 

where Mq is a scaling mass and the profile factor, Vmas, has the explicit 
expression: 

l^mas^S TfiOed^ ; (5) 

Jo 

and the integration is carried up to: 

H = - ; (6) 

^0 

which may be conceived as a scaled, truncation radius. 

Two limiting cases of Eq. (3) can be relevant to the aim of the current 
paper. The first one occurs as o; — > +00, which implies x provided /3 
and 7 remain finite, and: 

+ ; e<i ; 



accordingly, Eq. (3) reduces to: 



= i-!' (7) 



which is a double power-law with exponents equal to —7 and —f3. 

The second limiting case of Eq. (3) occurs as /3 — > +00, which implies 
X +00 provided a and 7 remain finite. If, in addition, ro —>■ +00 i.e. 
S — > 0, and: 



where and Cp are constants, then Eq. (3) reduces to: 

lim /(O = § ; (8) 



A similar result is also obtained as a — > provided x remains constant, which 
implies p ^ J and: 

accordingly, Eq. (3) reduces to: 

\im f{0^2-^r' ; (9) 

a— >U 

which is a power-law with exponent equal to —7. 

The logarithmic GPL density profile, deduced from Eq. (3), is: 

iog/ = -7ioge-xiog(i+n ; (10) 

and the maximum variation in slope occurs at log^ = log(r/ro) = i.e. 
r = Tq, which discloses the geometrical meaning of the scaling radius in GPL 
density profiles. In the logarithmic plane (0 log^ logf), Eq. (10) represents 
a curve where two asymptotes exist: 

iog/ = -7ioge ; e«i ; (iia) 
iog/ = -/3ioge ; e»i ; (lib) 

which meet at the origin. For further details, see GMV. 



2.2 SDH density profiles 

Simulated haloes are characterized by a "virial" parameter, either the virial 
mass, Myir, or the virial radius, Vyir, defined such that the mean density 
inside the virial radius is A^j^ times the mean matter universal density, ph — 
Pcrit^m, at that redshift: 

47r 3 

■^vir — "^^virPcrit^mTmr ' (1^) 

where pcrit — (SttG) is the critical density for closure. The critical over- 
density at virialization, A^j^, is motivated by the spherical collapse model: 
it is below two hundreds for an Einstein-de Sitter cosmology, and exceeds 
three hundreds for a flat AGDM cosmology where ~ 0.3, at ^ = (e.g.. 
Bullock et al. 2001). 
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Plotting SDH density profiles on a logarithmic plane, (Ologr; logip), 
necessarily implies use of dimensionless coordinates, defined as: 

r p 

?7 = — ; ipiv) = — ; Ph^ Pcrit^m ; (13) 

"^vir Ph 

where an upper and lower limit to the domain follow from the definition of 
virial radius and the occurrence of numerical artifacts, respectively. This is 
why (i) regions outside the virial radius are still falling in, and then macro- 
scopic kinetic energy has still to be converted into peculiar energy (e.g.. Cole 
& Lacey 1996; NFW), and (ii) two-body relaxation takes place in the central 
regions of simulated haloes, within about O.Olryij. (e.g.. Bullock et al. 2001; 
Fukushige & Makino 2001, 2003, 2004; Navarro et al. 2004). Accordingly, 
SDH density profiles within the range: 

-2 < logr; < ; (14) 

appear to be closely related to the virialized region. 



2.3 Comparison between GPL and SDH density pro- 
files 

A comparison between GPL and SDH density profiles necessarily implies that 

the truncation radius, and the mass enclosed within the truncation isopycnic 
surface, do coincide with the virial radius and the virial mass, R — ^vir si^nd 
M = Myir, respectively. Then the combination of Eqs. (4), (12), and (13), 
yields: 

Pfi ^mas 

^vir ) (1^) 

^0 

where the scaled virial radius, ^^ir, may be taken as a definition of the con- 
centration (CMV). 

Relating scaled GPL density profiles, expressed by Eqs. (2), to scaled SDH 
density profiles, expressed by Eqs. (13), yields: 

^^Urv ; / = -^ ; (17) 

Po 
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where is defined by Eq. (16). Accordingly, a generic, scaled GPL density 
profile, expressed by Eq. (3), takes the equivalent form: 

and the related, logarithmic GPL density profile, is deduced by use of Eq. (15). 
The result is: 

log ^ = log A^ir - log Umas + 3 log^„i^ 

-7 log Cvir - 7 log 7/ - X log [1 + iCvirVT] (19) 

which depends on three exponents, (a,/3, 7), and two scaling parameters, 
(^O)Po)- On the other hand, the scaled mass, i^mas, is defined by Eq. (5); 
the virial radius, Vyir, is known from the computer run; and the critical 
overdensity, A^jj^, is determined by the cosmological model. 

In the logarithmic plane, (0 log rj log ■0) , Eq. (19) represents a curve where 
two asymptotes exist: 

log 1p = log Ayir - log Umas + 3 log ^yir 

-7l0gC^;ir- - 7l0g?7 ; Cmr-?7<1 ; (20a) 

log 1p = log Ayir - log Umas + 3 log ^yir 

-/^logCmr- -/31og77 ; ^„i^77 > 1 ; (20b) 

which meet at the point (log 77, log ■0) = [log(ro/rOTr), log(po/P/i)]- For further 
details, see CMV. 

2.4 Analytical and geometrical properties of GPL den- 
sity profiles 

The change of coordinates: 

a; = log77 ; y = \og^jJ ; (21) 

makes the expressions of the curve and the asymptotes, defined by Eqs. (19), 
(20a), and (20b), translate into: 

y = yc - l{x - xc) - x'^og{l + expio[a{x - xc)]} ; (22) 
y = yc - lix - xc) ; x <^xc ; (23) 
y = yc - l3{x - xc) ; x » ; (24) 
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where, in general, exp^a; = , and expo; = e^, according to the standard 
notation. The point, C{xc, yc), represents the intersection of the asymptotes 
in the plane (Oxy), and the expressions of the coordinates are (CMV): 

xc^ - log ^yir = - log — = log — ; (25) 
g = log— ; (26) 

^mas Ph 

which allow the following expression of the scaled mass: 

log Vmas = log ^vir " 3xc " yc ■ (27) 

The vertical intercept of the curve, 6, and the asymptotes, and bp, may be 
determined through the combination of Eqs. (22), (23), (24), and (27). The 
result is: 

b = yc + l^c - X log [1 + expio(-aa;c)] ; (28) 
b^ = yc + ixc ; (29) 
b(3 = yc + f3xc ; (30) 

for further details, see CMV. 

It can be proved (CMV) that for arbitrarily chosen SPH density profile, 
cosmological model, and GPL density profile, there is a one-to-one correspon- 
dence between analytical parameters, (ro, po, a, (3, 7), and geometrical param- 
eters, {xc,yc,b,bi3,by), in the sense that either set univocally determines a 
GPL density profile. Owing to Eqs. (25), and (26), the set {xc,yc,Oi, 
may also be used to this respect. 



3 An analytical RFSM5 method 

Let {xi,yi) be coordinates of a generic point of a logarithmic SDH density 
profile, and [xi,y{xi)] be their counterparts related to a fitting, GPL density 
profile. Owing to Eq. (22), the corresponding, logarithmic absolute residual, 
is: 

/? - 7 

Ri = {yi - yc) + - xc) H log{l + exp^Q[a{xi - xc)]} ; (31) 

a 
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where the exphcit expression of the exponent, x, has been used, according 
to Eq. (1). 

With regard to a SDH density profile, a RFSM5 (Residual Function Sum 
Minimization within a 5-dimension hypcrspacc) method consists in the min- 
imization of the sum - over the range defined by Eq. (14) - of an assigned 
residual function (e.g., square, R^, or absolute value, \R\), within a selected 
domain of a 5-dimension hyperspace, (0 Xc Vc P j), or any equivalent 
choice of parameters. A numerical procedure has been used in CMV, where 
the sum of absolute values and squares of residuals was performed on 10^ 
nodes of a hypergrid within the above mentioned 5-dimension hyperspace. 
On the contrary, an analytical procedure shall be exploited here. 

In general, let v = F{ui,U2, ■■■,un) be a function which is defined in a 
domain of a A^-dimension hyperspace, {OuiU2-..uisf), and continuous therein. 
Then a necessary and sufficient condition for a point, P*{ul,U2, be- 
longing to the domain, to be an extremal mimimum, is that the coordinates 
of the point under discussion make (i) a solution of the system, where all the 
first partial derivatives equal zero: 

and (ii) the Hessian determinant positive: 

>0; 1<J<A^; 1<A;<A^; (33) 



^dujduk J 

in the more general case where at least one first partial derivative, {dF/ duj)p* , 
docs not exist, Eq. (32) makes a necessary condition for a point, P*{ul, 
u*j^), to be an extremum point of minimum. For further details see e.g., 
Smirnov (1969, Chap.V, §§22-26). 

Owing to the discontinuity of the first derivative of the absolute value, 
for null arguments, which makes a solution of the problem much more dif- 
ficult, our attention shall be limited to the sum of square residuals in the 
5-dimension hyperspace of interest: 

n 

F{xc,yc,a,^,^) = Y.Rl ; (34) 
1=1 
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where the i-th residual, Ri, is defined by Eq. (31). Accordingly, Eq. (32) 
reduces to: 

where Ui = xc, U2 = He, ""s = a, = [3, = 7. 

The combination of Eqs. (31), (34), and (35), where the first derivatives 
are written explicitly, after long but stimulating algebra yields three alterna- 
tive expressions of the exponent, 7, as: 

Vs{x,s)vs{y,t) -Vs{s,t)vs{x,y) 
7 = 71 xc, a) = -, . , , , r ; 36a) 

Vs{x, S)Vs{x, t) - Vs{s, t)Vs{x, X) 

7 = l2[xc, a) = ^ ; 36b 

Vs{x,s)vs{y,xt) - Vs{s,xt)vs{x,y) 
7 = l3[xc, a = 7 ^-7 7 rr-T r ' 

Vs{X, S)Vs[X, Xt) — Vs[S, Xt)Vs{X, X) 



where, in general: 



Vs{a, h) — ah — ah ; (37a) 
1 



c 



— ^ Q ; c — a,h,ah; (37b) 



and the remaining variables, Sj and ti, are defined as: 



Si = log(l + Wi) ; (38a) 
ti-TJ^ ; (38b) 
Wi = expio[Q;(xj - Xc)] ; (38c) 



where, in addition: 



dxc 

for a formal demonstration, see Appendix A 



= Wi{xi - Xc) In 10 ; (39a) 
= -WialnlO ; (39b) 
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In a 3-dimension space, (0 xc ol 7), Eqs. (36) represent three surfaces. 
Let {x}j a* 7*) be an intersection point of the above mentioned surfaces. The 
remaining parameters, /3* and y^, may be calculated using Eqs. (61) and (59), 
respectively. The resulting point in a 5-dimension space, (0 xc yc P 
let it be P* = (x^, y^, a*, 7*), satisfies the necessary condition for the 
occurrence of an extremal minimum, expressed by Eqs. (48). Then a key role 
is played by the intersections of the above mentioned, three surfaces. 

To this regard, the following algorithm is used. First two pairs of surfaces 
are considered, let they be (71,72) and (71,73), then the intersecting lines 
related to each pair are simultaneously determined. If the intersecting lines 
exhibit an intersection point, then the necessary condition for the occurrence 
of an extremal minimum is satisfied at that point. In the meantime, it 
can be numerically tested if the intersection point is, in fact, a minimum, 
thus avoiding the calculation of the Hessian determinant. On the contrary, 
nothing can be said about the existence of other extremal minima, and the 
coincidence of the point under consideration with an absolute minimum. 

Finding no intersection point does not necessarily mean absence of ex- 
tremal minima: it could imply that the algorithm used here is insufficient 
to this respect. If it is the case, a non extremal minimum is selected in the 
subdomain where the algorithm works. 

4 An application to SDHs on the scale of 
cluster of galaxies 

Using an analytical RFSM5 method, GPL density profiles arc fitted to a 
sample of 17 SDH density profiles, on the scale of cluster of galaxies within 
a fiat ACDM cosmology (Rasia et al. 2004), as already done in CMV. The 
values of the cosmological parameters used therein are: Qa = 0.7; Jl^ = 0.3; 
Q,}, — 0.03; h = 0.7; as = 0.9; where the symbols have their usual captions 
(e.g., Klypin et al. 2001; Bullock et al. 2001) and, in particular, the indices 
m and b denote all matter (dark + baryons) and baryons only, respectively. 
For a detailed discussion on the computer code, initial conditions, resolution 
issues, and the way of finding the halo centre, see Rasia et al. (2004) and 
further references therein. 

Simulations include both dark and baryonic matter, but only the former 
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is relevant to the aim of the current paper. Accordingly, the baryonic matter 
shall not be considered, and all the parameters shall be intended as related 
to the dark matter halo. 

The definition of the virialized region within each halo, via Eq. (12), needs 
the knowledge of the critical overdensity at virialization, A^ir. With regard 
to total (dark + baryonic) matter, it depends on the cosmological model 
(e.g.. Bullock et al. 2001) and, in the case under discussion, (A^j^)^ = 323 
at 2; = 0, where all the sample haloes may safely be thought of as viriahzed 
(Rasia et al. 2004). 

If only the dark matter is considered, then {Amr)cdm = A^ir = C{A.mr)m, 
where ( is the fraction of dark matter in each density perturbation, and 
averaging over the whole sample yields ( = 0.907 (Rasia et al. 2004). Ac- 
cordingly, the value: 

A^ir = 0.907 • 323 = 292.961 ; (40) 

can be used to an acceptable extent^. 

Owing to the random criterion used for the selection (Rasia et al. 2004), 
sample haloes are characterized by varying dynamical properties: at the 
present time, some are more relaxed, while others are dynamically perturbed. 
The surrounding environment can also be quite different: some clusters are 
more isolated, while others are interacting with the surrounding cosmic web. 
Accordingly, the related results shall be representative of an averaged cluster, 
in an averaged environment and dynamical configuration. 

4.1 SDH density profiles 

The main features of sample haloes at z = 0, are listed in Tab. 1. The mass, 
been taken from Rasia ct al. 2004, while the mass, Al',-,., has been 
deduced from Eq. (12). The apparent discrepancy between Myir and M^^ is 
owing to two different sources. First, a systematic contribution takes origin 
from the uncertainty on A„jr and, second, a random contribution arises from 

^The above value of the critical overdensity was deduced from an earlier, unpublished 
version of Rasia ct al. (2004). It is slightly different from A„,r = 0.903 x 323.7625 = 
292.3576 deduced from the current version, which appeared when all the calculations 
of this paper were performed. As the relative difference amounts to about 0.2%, the 
calculations were not repeated using the latter value of the critical overdensity. 
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427583 
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49359 
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28433 
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68262 
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36489 
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159330 
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71690 


71463 


15 
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107229 
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48250 


47971 
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58734 


1524 


36060 


36132 


17 


510.07 
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71937 


1628 


44170 


44045 



Table 1: Main features of sample haloes at 2; = 0. Column captions: 1 - case; 
2 - computer run; 3 - type {R - safely relaxed; M - safely a major merger 
occurring); 4 - number of dark matter particles within the virial radius; 5 - 
virial radius {h~^ kpc); 6 - virial mass (/i~^lO^°m0); 7 - virial mass deduced 
from Eq. (12). Both virial radii and virial masses are normalized to the 
dimensionless Hubble parameter at the current time, h. 
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the uncertainty on Vmr. An additional random contribution is related to 
averaging the fraction of dark matter over the whole sample, with regard to 
Eq.(12). 

The relative difference, |1 — M^^/Myj^l, is less then one percent in all cases 
except 9, where it is less than one and half percent. The virial mass can be 
evaluated, to a good extent, by use of Eq. (12), taking the virial radius from 
the results of simulations. It is worth mentioning that the RFSM5 method 
is independent of the value of the virial mass, while a change in virial radius 
makes SDH density profiles systematically shift along the horizontal axis in 
the logarithmic plane (0 log?] logip). 

The mean SDH density profile is obtained by averaging over the whole 
set the values related to each logarithmic radial bin, in the range of interest, 
expressed by Eq. (14). The value of the critical overdensity at viriahzation, 
A^jr; is fixed by Eq. (12), then a single free parameter remains: the virial 
radius, r^ir, which allows the calculation of the virial mass, M^ir. To this 
aim, first the virial radius is averaged over the whole sample, using the data 
listed in Tab. 1, and then the mean value of the virial radius is inserted into 
Eq. (12). The result is: 

r^ir = 1866 h-^ kpc ; M^ir = 66330 h'^ lO^^m© ; (41) 

and the apphcation of a RFSM5 method yields best fitting GPL density 
profiles, with radius and mass equal to f^r and M^ir, respectively. 

A mean virial radius has been preferred in place of a mean concentration 
(Rasia et al. 2004) for the following reasons. First, virial radii are indepen- 
dent of GPL density profiles, contrary to concentrations, or velocity profiles, 
which should be calculated for any choice of the fitting profile. Second, the 
range in virial radius, 1407 < r^ir/ (h^^kpc) < 2470, corresponds to relative 
errors of about 25% and 33%, respectively, with regard to a mean value, 
r^r = 1866 h'^ kpc. On the other hand, the range in concentration (cal- 
culated for NEW density profiles), 5 < ^^j^ < 10, corresponds to relative 
errors of about 32% and 37%, respectively, with regard to a mean value, 
^vir — 7.2976. Then the average of the virial radius should be preferred to 
this respect. 

Having in our hands a sample of SDH density profiles, listed in Tab. 1, 
together with the averaged virial radius, and related virial mass, expressed 
by Eqs. (41), we are left with the search of a best fitting GPL density profile. 
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which minimizes the sum of square residuals. To this aim, the analytical 
RFSM5 method formulated in section 3, shall be used. 

4.2 Results 

As outlined in section 3, an analytical RFSM5 method has been applied to 
each sample halo, hsted in Tab. 1, and to the related, mean SDH density 
profile, which has been defined above. The values of the exponents, a, x, /3, 
7, the scaled radius, ^^r, the scaling radius, vq, the scaling density, po, and 
the sum of absolute value residuals, 1-^*1) ^^'^ square residuals, J^Rh 
the extremal or fiducial minimum, are listed in Tabs. 2 and 3, respectively, 
for both analytical (current paper, upper lines) and numerical (CMV, lower 
lines) RFSM5 methods. The following conclusions are deduced. 

(i) A necessary condition for the occurrence of an extremal point, expressed 
by Eqs. (36), is satisfied for 8 SDHs and is not for the remaining 9 
together with ADP. In the former alternative, an extremal minimum 
point may safely exist. In the latter alternative, the existence of an 
extremal minimum point cannot be established, but only a non ex- 
tremal minimum may be determined. This last is related to the GPL 
density profile where the distance between intersecting fines of two dis- 
tinct pairs of surfaces, defined by Eqs. (36), attains the lowest value. 
The above mentioned, non extremal minimum, may be located both 
within (IMP) and on the boundary (BMP) of the subdomain where 
the RFSM5 method, used in the current paper, works. 

(ii) The sum of square residuals is systematically lower using analytical 
(current paper) instead of numerical (CMV) RFSM5 methods, provided 
an extremal minimum point occurs. Accordingly, EMPs may safely be 
thought of as absolute minima. With regard to non extremal minimum 
points, the above result maintains in 4 cases (including ADP) while the 
contrary holds for the remaining 5 cases. The trend is not necessarily 
the same for the sum of absolute value residuals, as expected. 

(ni) The exponents, a, /3, 7, range within 0.4 < a < 19.8 {a — 0.2900 
for ADP); 1.9 < p < 17.6; 0.6 < 7 < 1.8; and 0.4 < a < 0.9; 
2.5 < f3 < 3.9; 0.7 < 7 < 1.3; using analytical and numerical RFSM5 
methods, respectively. With regard to SDH density profiles where 
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Table 2: Exponents of GPL density pr|a|iles which minimize the sum of square 
residuals, using (i) an analytical RFSM5 method (current paper, upper lines) , 
and (ii) a numerical RFSM5 method (CMV, lower lines), with regard to 
17 sample haloes listed in Tab. 1, together with the averaged density profile 
(ADP). Type captions: I - safely relaxed; M - safely a major merger occurring. 
Class captions: EMP - extremal minimum point; IMP - minimum point 
occurring within the subdomain interested by the application of the method. 













Tn 


























1 


T 

J. 


s 

(J 


5258 


hj—z 




2724 


E-l-4 


1 

X 


4808 

t:0 WO 


E- 


_7 


2 


755o 


tl — 


1 

X 


XJ 


0097 


Hi- 


-3 

XJ 


IMP 








DoUo 


rj+U 


7 


6843 




4 


5733 


E- 

XJ 


XJ 


4 


oi on 


bj— 


1 

X 


1 

X 


z84o 


H- 


-2 




2 


M 


1 


(Jzoo 




(J 


2172 

^ X f ^ 


E+2 


5 


1755 

X f \J\J 


E- 


\J 


fi 


4089 


hi— 


1 


3 


o7d9 


H- 


-2 


BMP 






7 


i4oU 




4 


6086 


E-l-2 


1 

X 




E- 

XL/ 


-4 


6 




hj— 


1 

X 


K 

XJ 


nno /I 
1)984 


TP 

t- 


-2 






M 

IV J. 


2 


uuyy 


77 1 1 


1 

± 


1 038 


E-l-3 


2 


0053 


E- 


-4 


6 

\J 


9707 


tl — 


1 

X 


7 


rio AC\ 

5840 


TP 

Hi- 


-2 


IMP 

XIV J.± 






A 


4ooo 


77 1 r\ 

b+0 


4 


9668 


E-l-2 


1 

X 


1 755 


E- 


-4 


8 

o 


z81o 




1 

X 


XJ 


9798 


TP 

Hi- 


-2 




4 


T 

X 


2 


0090 


TT^ 1 1 


1 

X 


4333 


E-l-2 


2 


002Q 


E- 

XJ 


-4 


7 


ood7 


bj— 


1 

X 


XJ 


1 TO /I 

17o4 


H- 


-2 


IMP 






A 


1 1 1 c; 

iilo 


rj+U 


« 

o 


1548 


E-l-2 


2 


13Q1 


E- 

XJ 


-4 


7 


Mil 


tl— 


1 

X 


4 


/I OO /I 

48o4 


H- 


-2 






I 


7 


UOOD 


77 1 

b+0 


4 


0324 


E+2 


3 


0094 


E- 




3 


9790 




1 


1 


3080 


TP 

Hi- 


-2 


EMP 






4 


7ooo 


77 1 r\ 




9426 


E-l-2 


^ 


0721 


E- 

XJ 


-4 

t: 


4 

t: 


5994 


TP 


1 

X 


1 

X 


5935 


TP 

H- 


-2 






T 

J. 


2 


U44i 


1 1 


1 

X 


2999 


E-l-2 


1 

X 


7204 


E- 

XJ 


-4 


3 


(009 


ill- 


1 

X 


1 

X 


O^^ /I 7 

zo4 1 


Hi- 


-2 


EMP 

XJ1VJ.J. 






1 


DoUi 


17 1 1 


1 


6103 

WX \J(J 


E+2 


7 

f 


3632 


E- 


-4 




zzz4 


bj— 


1 


3 

xj 


o055 


H- 


-2 




7 

f 


M 


2 


1719 


TT^ 1 1 


1 


4451 


E+2 


3 

KJ 


3786 


E- 


-4 


7 

f 


9368 


Hi- 


1 


4 


8776 


H- 


-2 


IMP 






4 


r\r\r"y 

zzbl 


77 1 r\ 


7 


4256 


E-l-2 


1 

X 


7070 


E- 

XJ 


-4 


6 


1250 


TP 


1 

X 


2 


n 1 o T 

9137 


TP 

H- 


-2 




s 

O 


M 




A'T A1 

4741 


77 1 


6 


4459 


E-l-3 


1 

X 


9954 


E- 

XJ 


-4 

t: 


2 


AO CO 

4863 


TP 

ill- 


1 

X 


4 


conn 

5o9U 


TP 

Hi- 


-3 


IMP 






4 


Ibo/ 


77 1 n 


s 

O 


4645 


E-l-2 


1 

X 


31 90 


E- 

XJ 


-4 


2 


z44d 


TP 


1 

X 


4 


9oD/ 


Tp 

H- 


-3 






M 


Q 


4zii 


17 1 A 


2 


7234 


E+2 


2 


0100 


E- 


-4 


2 


6901 


Hi — 


1 


XJ 


/4dz 


H- 


-3 


EMP 






(J 


Dzo4 


17 1 A 


4 


5626 


E+2 


2 


2399 


E- 


-4 


2 


o09z 


Hi- 


1 


fi 


1 non 
iU89 


H- 


-3 




10 






3231 


77 1 -1 


7 


2696 


E+2 


3 


8804 


E- 


-4 


4 


6024 


TP 

t — 


1 


2 


4io9 


TP 

Hi- 


-2 


IMP 








bo (4: 


77 1 n 




6395 


E-l-2 


2 


51 32 


E- 

XJ 


-4 


c; 


/6o 


TP 


1 

X 


2 


nTo 1 
U Mi 


TP 
t- 


-2 




1 1 

X X 






no 1 


77 1 r\ 


2 


8983 


E-l-2 


3 


4633 


E- 

XJ 




3 


c)oyo 


t — 


1 

X 


9 


4Udz 


Tp 

H- 


-3 


EMP 

XJ1VJ.J. 






3 


/ lo4 


17 1 n 




4100 


E+2 


fi 


7645 


E- 




4 


5841 


Hi — 


1 


1 


/ D/4 


H- 


-2 




12 




1 




1 1 

rj+i 


1 


6432 


E+2 


1 


2941 


E- 


-4 


4 


809o 


H— 


1 


2 


OO yl /I 

8o44 


H- 


-2 


EMP 








i2Ui 


77 1 r\ 

b+0 


3 


6319 


E+2 


2 


8856 


E- 


-4 


8 


A A AO 

1443 


TP 

t — 


1 




r\ AOT 

0487 


TP 
t- 


-2 




X (_l 




8 

(J 


n 1 /I 
Do 14 


17^ 1 n 
hv+O 


2 


4964 


E+2 




2008 


E- 

XJ 




6 


lolo 


H— 


1 

X 


3 


on V/^ 
zy 


TP 

t- 


_2 


EMP 

XJ1VJ.± 








8234 


E+0 


3 


2012 

X^ 


E-l-2 


2 


5718 

1 XO 


E- 

XJ 


-4 




9097 


E- 


1 

X 


3 


6570 


E- 


-2 




14 


I 


1 


3437 


E+0 


2 


0338 


E+3 


3 


9969 


E- 


-5 


3 


8299 


E- 


1 


1 


1169 


E- 


-2 


EMP 






5 


3088 


E+0 


5 


1477 


E+2 


7 


7666 


E- 


-5 


4 


1735 


E- 


1 


1 


2273 


E- 


-2 




15 




8 


9712 


E+0 


2 


6673 


E+2 


4 


6219 


E- 


-5 


4 


9331 


E- 


1 


2 


0331 


E- 


-2 


EMP 






7 


6913 


E+0 


3 


nil 


E+2 


1 


9963 


E- 


-4 


6 


0948 


E- 


1 


2 


8821 


E- 


-2 




16 




5 


2999 


E+0 


4 


1079 


E+2 


3 


8564 


E- 


-5 


6 


4586 


E- 


1 


2 


9673 


E- 


-2 


IMP 






3 


7154 


E+0 


5 


8599 


E+2 


1 


1226 


E- 


-4 


4 


4819 


E- 


1 


1 


5259 


E- 


-2 




17 




1 


4047 


E+1 


1 


6556 


E+2 


7 


2820 


E- 


-5 


7 


4326 


E- 


1 


7 


3645 


E- 


-2 


IMP 






4 


6774 


E+0 


4 


9723 


E+2 


1 


5144 


E- 


-4 


7 


1010 


E- 


1 


5 


7718 


E- 


-2 




ADP 




2 


2073 


E+0 


1 


2077 


E+3 


5 


9044 


E- 


-6 


1 


1027 


E- 


1 


1 


1783 


E- 


-3 


IMP 






3 


6308 


E+0 


7 


3422 


E+2 


1 


0238 


E- 


-4 


1 


3262 


E- 


1 


1 


4983 


E- 


-3 





Table 3: Scale parameters of GPL dgigsity profiles which minimize the sum 
of square residuals, using (i) an analytical RFSM5 method (current paper, 
upper lines), and (ii) a numerical RFSM5 method (CMV, lower lines), with 
regard to 17 sample haloes listed in Tab. 1, together with the averaged density 
profile (ADP). The scaling radius, ro, and scaling density, po, are expressed 
in kpc and lO^°m0/kpc^, respectively. Other captions as in Tab. 2 



an EMP occurs, the above ranges in the former alternative reduce to 
0.5 < a < 19.8; 2.0 < (3 < 4.4 (/3 = 4.4590 for ADP); 0.8 < 7 < 1.8; 
which is, only for P, substantially closer to their counterparts related 
to the latter alternative. Accordingly, analytical RFSM5 methods yield 
GPL density profiles where the range in the exponents is larger than in 
their counterparts related to numerical RSFM5 methods. It provides 
additional support to the occurrence of a certain degree of degeneracy 
(e.g., Klypin et al. 2001; CMV). 

(iv) The scaled radius, ^yir, the scaling radius, Tq, and the scaling den- 
sity, po, range within 0.08 < ^^ir < 33.3; 72 < ro/kpc < 32724; 
0.01 < lOVo/(lO^°m0/kpc^) < 3.9; and 3.6 < ^^ir < 16.5; 161 < 
ro/kpc < 847; 0.4 < lOVo/(lO^°m0/kpc^) < 7.4 ; using analytical 
and numerical RFSM5 methods, respectively. With regard to SDH 
density profiles where an EMP occurs, the above ranges in the for- 
mer alternative reduce to 1.3 < ^vir < 20.5; 129 < ro/kpc < 2034; 
0.3 < 10Vo/(10^°mo/kpc^) < 2.1; which is much closer to their coun- 
terparts related to the latter alternative. 

The occurrence of high values of the exponents, a and related to (a) 
lack of an EMP, and (b) lower value of the sum of square residuals using 
analytical instead of numerical RFSM5 methods, as in cases 1 and 10 of 
Tab. 2, suggests further investigation on the best fitting, limiting GPL density 
profiles, analysed in subsection 2.1. 

According to Eq. (7), GPL density profiles reduce to a double power-law 
(DPL) with exponents equal to —7 and — /3, in the limit a — > 00, x — > 0, 
i.e. 7 and f3 finite. Then we are left with two exponents instead of three, 
and a RFSM4 method has to be used, within a 4-dimension hyperspace 
(0 xc He P 7)) in dealing with each sample halo, listed in Tab. 1, together 
with the related ADP, which has been defined above. The values of the 
exponents, /3 and 7, the scaled radius, ^^r, the scaling radius, ro, the scaling 
density, po, the sum of absolute value residuals, and the sum of square 

residuals, J2R^, are listed in Tab. 4. By comparison with Tabs. 2 and 3, the 
following conclusions are deduced. 

(v) The sum of square residuals is systematically larger with respect to 
analytical RFSM5 methods, provided an EMP occurs, as expected. 
Concerning non extremal minimum points, the above result maintains 
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Table 4: Parameters of double power-law (DPL) density profiles which min- 
imize the sum of square residuals, using an analytical RFSM4 method, with 
regard to 17 sample haloes listed in Tab. 1 together with the averaged density 
profile (ADP). Other captions as in Tabs. 2-3. Cases, types, and classes have 
been omitted to save space. 
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in 6 cases (including ADP), while the contrary holds for the remaining 
4 cases. The sum of square residuals is larger, with respect to numerical 
RSFM5 methods, in 11 cases (including ADP), while the contrary holds 
for the remaining 8 cases. The trend is not necessarily the same for the 
sum of absolute value residuals, as expected. 

(vi) The exponents, (3 and 7, range within 1.9 < /3 < 3.5 and 1.1 < 7 < 
1.9, respectively, which are comparable to their counteparts related to 
both analytical (restricted to occurring EMPs) and numerical RSFM5 
methods. 

(vii) The scaled radius, ^mr-, the scaling radius, ro, and the scaling density, 
Po, range within 1.9 < < 33.3; 85 < ro/kpc < 1471; and 0.01 < 
lO'^Po/ (lO^°m0/kpc^) < 3.9; respectively, which are comparable to their 
counterparts related to both analytical (restricted to occurring EMPs) 
and numerical RFSM5 methods. 

Then DPL density profiles may be considered as a rough, but viable approx- 
imation in fitting SDH density profiles. 

According to Eqs. (8) and (9), GPL density profiles reduce to a simple 
power-law (SPL) with slope equal to —7, in the limit (3 — > +00 where a, 7, 
and (1 -|- ^°')^ remain finite, and a — >■ where (3 ^ 'j and (1 -|- ^°')^ — >■ 2^. In 
any case, we are left with one exponent instead of three and one coordinate 
(the vertical intercept of SPL density profiles in the logarithmic plane, b^) 
instead of two. Then a RFSM2 method i.e. a standard least-squares fit has 
to be used, within a plane (0 7), in dealing with each sample halo, listed 
in Tab. 1, together with the related ADP, which has been defined above. 
The values of the exponent, 7, the intercept, 6^, the sum of absolute value 
residuals, J2\Ri\: a-nd the sum of square residuals, J2 Rh listed in Tab. 5. 
By comparison with Tabs. 2 and 3, the following conclusions are deduced. 

(viii) The sum of square residuals is systematically larger with respect to 
both analytical and numerical RFSM5 methods. The same holds also 
for the sum of absolute value residuals. 

(ix) The exponent, 7, and the intercept, bj, range within 1.7 < 7 < 2.1 
and 2.0 < b^ < 2.2, which is larger and lower, respectively, with regard 
to their counteparts related to both analytical and numerical RSFM5 
methods, as expected. 
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Table 5: Parameters of single power-law (SPL) density profiles which mini- 
mize the sum of square residuals, using an analytical RFSM2 method i.e. a 
standard least squares fit, with regard to 17 sample haloes listed in Tab. 1 to- 
gether with the averaged density profile (ADP). Other captions as in Tabs. 2 
and 3. 
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Then SPL density profiles cannot be considered as an acceptable approxima- 
tion to GPL density profiles, in fitting SDH density profiles. 

To get further insight, the sum of square residuals is hsted in Tab. 6 for 
all the four alternatives which have been exploited in the current section. 
The values of the exponents, a, X) /^i 7) the scaled radius, C,vir, the scaling 
radius, Tq, the scaling density, po, and the sum of absohitc vahie residuals, 
X) l-Ril, and square residuals, J2Rh the extremal or fiducial minimum, are 
listed in Tabs. 7 and 8, respectively, with regard to the best fit, among the 
ones considered above, for each sample halo. 

The values of some analytical and geometrical parameters, tjadp, related 
to the best fitting GPL density profile to the mean SDH density profile, 
are listed in Tab. 9 together with their counterparts, rj, averaged over the 
best fitting GPL density profiles to the whole halo sample, using the data 
of Tab. 8. Also listed in Tab. 9 are the related standard deviations, (Tg ^, the 
standard deviations from the mean, o-g fj, and the standard deviations from 
the standard deviation from the mean, ag ji, which are expressed as (e.g., 
Oliva & Terrasi 1976, Chap. V, §5.6.3): 



1 " 



' STj 
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sri 
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nl/2 



1 1 " 



n n 



' srj 



— j- Y^ivi - v) 
^ i=i 

n 
i=l 



1/2 



(42) 

(43) 

(44) 
(45) 



where r] = a, p, 7, ^^ir, Vmas, log^^j^, logUmas, yc, b, hp, by, n = 17 in all 
cases but n = 14 if 77 = a, as a — +00, related to DPL density profiles, 
cannot be considered in averaging; and, owing to Eq. (25), log^^j^ = —xc- 
The following conclusions are deduced from Tab. 9. 



X 



Values of parameters, rjADp, related to the best fitting GPL density 
profile to the mean SDH density profile, are different from their coun- 
terparts averaged over the best fitting GPL density profiles to the whole 
halo sample, listed in Tabs. 7 and 8. The discrepancy changes for differ- 
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Table 6: The minimized sum of square residuals which provides the best 
fit to 17 sample haloes listed in Tab. 1 together with the averaged density 
profile (ADP), with regard to the following alternatives: numerical RFSM5 
methods involving GPL density profiles (GPN); analytical RFSM5 methods 
involving GPL density profiles (GPA); analytical RFSM4 methods involving 
DPL density profiles (DPA); analytical RFSM2 methods involving SPL den- 
sity profiles i.e. standard least squares method (SPA). Other captions as in 
Tabs. 2 and 3. 
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Table 7: Exponents of GPL or DPL density profiles which minimize the 
sum of square residuals, using both analytical (current paper) and numerical 
(CMV) RFSM5 or RFSM4 methods, with regard to 17 sample haloes hsted 
in Tab. 1, together with the averaged density profile (ADP). Method cap- 
tions: GPA - analytical RFSM5; GPN - numerical RFSM5; DPA - analytical 
RFSM4 where GPL density profiles are reduced to their DPL counterparts 
i.e. a — > +00 and x ^ 0- Other captions as in Tab. 2. 
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Table 8: Scale parameters of GPL or DPL density profiles which minimize the 
sum of square residuals, using both analytical (current paper) and numerical 
(CMV) RSFM5 or RSFM4 methods, with regard to 17 sample haloes listed 
in Tab. 1, together with the averaged density profile (ADP). Method cap- 
tions: GPA - analytical RFSM5; GPN - numerical RFSM5; DPA - analytical 
RFSM4 where GPL density profiles are reduced to their DPL counterparts 
i.e. a —>■ +00 and % — > 0. Other captions as in Tab. 2. Cases and types have 
been omitted to save space. 
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Table 9: Comparison between parameters, tjadp, related to the best fitting 
density profile to the averaged SDH density profile (ADP), with regard to 
analytical and numerical RFSM5 or RFSM4 methods, respectively, and their 
counterparts, rj, averaged over the best fitting density profiles to the whole 
halo sample, hsted in Tabs. 7 and 8. Also listed are the related standard de- 
viations, (Ts ^, the standard deviations from the mean, cTg ^, and the standard 
deviations from the standard deviation from the mean, ag ji- With regard to 
the exponent, a, three DPL density profiles have been excluded as a — +oo 
in that limiting situation. To save space, the notations, = log^„i,. and 
= log Uraas, have been used. It is worth remembering that log^^j^ = —^c, 
according to Eq. (25). 
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ent parameters, from less than =1=0-^ fj to more than ^Sus fj. The above 
result is expected from the theory of the errors. 

(xi) The exponents of best fitting, GPL density profiles related to ana- 

lytical RSFM5 methods, are distant from their NFW counterparts, 
(a,/3,7) = (1,3,1), conform to [Nint(a), Nint(/5), Nint(7)] = (0,4,2), 
where Nint denotes the nearest integer. The difference increases from 
about one half for 7 to about one and half for f3. On the other hand, the 
contrary holds with regard to GPL density profiles related to numerical 
RFSM5 methods, namely [Nint(a), Nint(/?), Nint(7)] = (1,3,1), where 
the difference increases from about one hundredth for 7 to about one 
half for a (CMV). 

(xii) The rsm error of the logarithm of the scaled radius, ^mr, is Cs iog{„i,. ~ 
0.60, to be compared with ag log^^i,. — 0.18 deduced from richer samples 
where (q;,/3,7) = (1,3,1); M^i/=(0.5-1.0) xlO^/i-^mo; 11 < n < 14; 
and n is an integer (Bullock et al. 2001). 

4.3 Discussion 

The application of an analytical RSFM5 method makes a necessary condition 
for the occurrence of an extremal minimum, Eqs. (36), which is satisfied in 
8 sample haloes on a total of 17, listed in Tab. 1. The related points in the 
5-dimension parameter hyperspace, P* = (a;^, |/^, a*, 7*), may safely be 
thought of as absolute extremal minima, by comparison with fiducial minima 
calculated using a numerical RSFM5 method (CMV). On the other hand, the 
existence of other extremal minimum points cannot be excluded in absence 
of a rigorous mathematical proof. In addition, the lack of extremal minima 
in the remaining 9 sample haloes plus ADP, could be an artifact related to 
the algorithm used for determining the intersection of the three surfaces, 
7 = '-f{xc,C(), expressed by Eqs. (36). An improved algorithm would suceed 
in localizing extremal minimum points in a larger number of sample haloes. 

The occurrence of an extremal minimum, or the algorithm power in find- 
ing it, seems to be correlated with halo morphology. The relative frequency 
of the above mentioned event is 7/12?»0.583 for safely relaxed haloes, and 
1/5=0.2 for safely undergoing a major merger haloes. 

The presence of a (safely absolute) extremal point of minimum, makes a 
firm criterion for deciding which, among the family of GPL density profiles. 
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best fits an assigned SDH density profile. In absence of extremal minima, 
numerical and analytical RFSM5 or RFSM4 methods, the last using DPL 
instead of GPL density profiles, seem to be competitive. With regard to 
sample haloes where no extremal minimum point has been found, the sum of 
square residuals is minimized using a numerical RFSM5 method in 3 cases, 
an analytical RFSM5 method in 2 cases, and an analytical RFSM4 method 
in 3 cases, as can be seen in Tab. 6. 

The related parameters, including cases where an extremal point of min- 
imum occurs, together with ADP, are shown in Tabs. 7 and 8. The ex- 
ponents, a, /?, 7, range within 0.42 < a < +00 (a = 0.29 for ADP); 
< X < 13 (x = 14.846 for ADP); 2.0 < (3 < 17.6; 0.6 < 7 < 1.9, 
respectively. The scaled radius, ^yir, the scaling radius, ro, and the scal- 
ing density, po, range within 0.08 < ^^ir < 20.5; 129 < ro/kpc < 32724; 
0.001 < lOVo/(lO^°m0/kpc^) < 2.6; respectively. The range in inner slope, 
7, is consistent with both core-shaped (7 < 1) and cuspy (7 > 1) density 
profiles The range in outer slope, f], is consistent with both infinite {j3 < 3) 
and finite (/3 > 3) total (with no truncation) mass. 

With regard to the exponent, 7, it is worth mentioning that, under reason- 
able boundary conditions such as force-free halo centre and vanishing density 
at infinite distance, the inequality 7 < 1 holds (Miicket & Hocft 2003). On 
the other hand, simple analytic considerations suggest that 7 < 2 (Williams 
et al. 2004). Finally, the validity of Jeans equation implies 1 < 7 < 3 
(Hansen 2004), where shallower slopes can arise, provided the effects of bary- 
onic matter are considered (El-Zant et al. 2004; Hansen 2004). The current 
result, 0.6 < 7 < 1.9, is consistent, or marginally consistent, with the above 
mentioned constraints. 

On the other hand, values of asymptotic inner slopes of fitting density pro- 
files determined in the current paper, are consistent with their counterparts 
deduced from recent high-resolution simulations using a three-parameter fit 
involving scaling radius, scaling density, and asymptotic inner slope (6 sam- 
ple objects, Dicmand et al. 2004) or a two-parameter fit involving scaling 
radius and asymptotic inner slope (16 sample objects. Reed et al. 2004). 
The related parameters are listed in Tab. 10, which shows agreement between 
different approaches, within the fiducial range, 7^3(7s ^. A marginal discrep- 
ancy between results from CMV and Reed et al. (2004) can be ehminated, by 

^Strictly speaking, a central cusp occurs for 7 > 0. 
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c n 7 as -y (Ts ^ A7 A+7 A^7 7 7+ 

T 17 1.3032 0.34209 0.082969 0.70036 0.43350 0.56693 1.0543 1.5521 

E 8 1.3218 0.28798 0.10181 0.45124 0.38076 0.41600 1.0164 1.6273 

C 17 1.0018 0.13893 0.033696 0.29868 0.19820 0.24844 0.90071 1.1029 

D 6 1.1617 0.13732 0.056060 0.24167 0.25833 0.25000 0.99352 1.3299 

R 16 1.2875 0.23910 0.059774 0.28750 0.41250 0.35000 1.1082 1.4668 

R 13 1.2692 0.24285 0.067353 0.26923 0.43077 0.35000 1.0671 1.4713 



Table 10: Comparison between statistical parameters related to the asymp- 
totic inner slope, 7, deduced from different samples using different fits. Col- 
umn captions: 1 - case: T - current paper, all sample haloes; E - current 
paper, only sample haloes where an extremum point of minimum occurs; C - 
Caimmi et al. (2004); D - Diemand et al. (2004); R - Reed et al. (2004); 2 - 
n: total number of sample objects; 3-7: arithmetic mean; 4 - CTg ^: standard 
deviation; 5 - -: standard deviation from the mean; 6 - A~7: maximum 
negative deviation from the mean; 7 - A+7: maximum positive deviation 
from the mean; 8 - A^7: mean maximum deviation from the mean; 9 - 
7~ = 7 — 3(Ts Ty- lower limit assigned to the mean; 10 - 7"^ = ^ + 3as ^: upper 
limit assigned to the mean. Sample haloes represent clusters or groups, with 
the exception of three objects belonging to the richer R sample (excluded in 
the poorer R sample), which represent SDHs embedding the Milky Way and 
two dwarf galaxies, respectively. 
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removing from the latter sample three objects, which represent SDHs embed- 
ding the Milky Way and two dwarf galaxies, respectively, instead of clusters 
or groups. On the other hand, current results make an improvement with 
respect to CMV. Mean values and related scatters, 7 ^ Scr^ ;y, corresponding 
to each case of Tab. 10, are represented in Fig. 1. 

No evident correlation is found between SDH dynamical state (relaxed or 
merging) and asymptotic inner slope of the logarithmic density profile, —7, or 
(for SDH comparable virial masses) scaled radius, ^^^r, contrary to previous 
results (Ascasibar et al. 2004) related to a sample of 19 high-resolution SDHs 
on the scale of both clusters of galaxies (13 objects) and galaxies (6 objects), 
with regard to NFW density profiles. An investigation on richer samples 
could provide more information to this respect. 

The large range in parameters shown in Tabs. 7 and 8, provides further 
evidence to the occurrence of a certain degree of degeneracy in fitting GPL to 
SDH density profiles (Klypin et al., 2001; CMV). For instance, with regard 
to case 14 listed in Tabs. 2 and 3, using analytical and numerical RFSM5 
methods, the sum of square residuals changes by less than 10%, while the 
remaining parameters change by more than 40%. 

As different methods are involved, it seems difficult that the occurrence 
of degeneracy is intrinsic to the procedure of minimization. A remaining 
possibility could be the lack of resolution in the central regions, which makes 
the lower limit of the domain, defined by Eq. (14). A more extended (on the 
left) range, owing to the next generation of high-resolution simulations, could 
reduce or eliminate the degeneracy. The fit has necessarily to be restricted to 
the virialized region of the halo, which leaves the upper limit of the domain 
unchanged. 

The occurrence of degeracy could also bias the procedure used for mini- 
mizing the sum of any function (in particular, the square) of residuals. For 
instance, the apphcation of a numerical RFSM5 method to the halo sample 
listed in Tab. 1 yields an inner slope, 7^1 (CMV), different from a value 
7 fa 1.5 currently found for dark matter haloes on the scale of cluster of galax- 
ies (e.g., Fukushige & Makino 2003; Hiotelis 2003), but closer to 7 f« 1.2-1.3 
from more recent attempts (Diemand et al. 2004; Reed et al. 2004). On the 
other hand, the apphcation of an analytical RFSM5 method to the above 
mentioned halo sample yields an inner slope, 7 ^ 1.5 (current paper). 

The scaled radius, = Tmr/f^o, is defined in the current paper as re- 
lated to the point where the maximum change in logarithmic slope occurs 
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Figure 1: Mean values and related scatters, 7=1=30-5 -, corresponding to cases 
listed in Tab. 10, positioned and labeled as therein. 
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in the corresponding GPL density profile (CMV). On the other hand, the 
concentration, c = rmr/'^'-2-i is usually defined as related to the point where 
the logarithmic slope equals —2 (e.g., Klypin et al. 2001; Bullock et al. 2001; 
Hiotelis 2003). The former definition seems to be more general, as it provides 
the scaled radius meaningful also in early times, where the slope of a GPL 
density profile may be smaller (in absolute value) than 2, while the contrary 
holds for the concentration (Hiotelis 2003). 

The GPL density profile which best fits to the averaged SDH density pro- 
file, with regard to analytical RFSM5 methods, is characterized by exponents, 
(a,/3,7), satisfying [Nint(a), Nint(/3), Nint(7)] =(0,4,2), which is different 
from NEW density profiles currently used in literature, (a,/?, 7) = (1,3, 1). 
The comparison with values averaged over the whole halo sample, discloses 
that fluctuations of the exponents around their NFW counterparts may safely 
be excluded. Accordingly, NFW density profiles cannot be conceived as in- 
trinsic to dark matter haloes of the current sample. 

On the other hand, following e.g.. Bullock et al. (2001), NFW density 
profiles (or alternative functional forms) may be considered as a convenient 
way to parametrize SDH density profiles, without implying that it necessarily 
provides the best possible fit. This is why the scaled radius, ^^r, and the 
scaled mass, Vmas^ can be interpreted as general structure parameters, not 
necessarily restricted to a specific density profile. In particular, any spread 
in ^vir and Vmas can be attributed to a real (larger) scatter in a "physical" 
scaling radius, defined by e.g., the radius where the change in slope of the 
logarithmic density profile attains its maximum value (CMV), rather than to 
inaccuracies in an assumed, "universal" density profile. For further details 
see e.g.. Bullock et al. (2001). 

Additional support to the above considerations is provided by the value 
calculated for the standard deviation of the decimal logarithm of the scaled 
radius, (Jsiog^^i^ ~ 0.59, which is 3.3-3.4 times larger than (J^iog^^j^ ~ 0.18 
deduced from a statistical sample of about five thousands simulated haloes, 
within mass bins equal to (0.5-1.0) xlO^/i^^m©, where 11 < n < 14 and 
n is an integer, with regard to NFW density profiles (Bullock et al. 2001). 
A comparable scatter, (Jsiog^^^^ ~ 0.17, has been obtained using numerical 
RSFM5 methods (CMV). The discrepancy with the results of the current 
paper, is probably owing to a bias in the selection of the subdomain where 
numerical RFSM5 methods work. 

To get further insight, the histogram of the decimal logarithm of the scaled 
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radius, log ^mr, is represented in Fig. 2, with regard to numerical (CMV) and 
best fitting (current paper, Tab. 8) RSFM5 methods. The former case (top 
panel) shows a clear cut at log^^jr — 0.55, which could be a biasing effect, as 
mentioned above, but the dispersion is low. The latter case (bottom panel) 
shows no sign of cut, but the dispersion is large. It is reduced only slightly by 
restricting the calculations to a homogeneous subsample, where only EMPs 
are considered, and even if the maximum deviation from the arithmetic mean 
is excluded therein. If a large dispersion still maintains for richer samples 
{N S> 17), in dealing with analytical RSFM5 methods, a low dispersion found 
in A^-body simulations (e.g.. Bullock et al. 2001) seems to be an artefact, 
related to the assumption of NFW (or any equivalent choice) density profile. 

On the other hand, a large spread in parameters of GPL density profiles 
where an extremum point of minimum occurs, may be considered as fiducial 
at most. This is in absence of a rigorous proof on the existence of an absolute 
point of minimum. If, in fact, absolute minima do not coincide with the 
extremum points of minimum which are found in the current attempt, then 
the spread in parameters could be reduced. 

The results of the current paper confirm a certain degree of degeneracy 
in fittig GPL to SDH density profiles, as pointed out by Klypin et al. (2001). 
For instance, cases 4 and 14 of Tabs. 2 and 3 exhibit - on the related upper 
and lower lines - values of the sum of square residuals and absolute value 
residuals, which are close each to the other. On the contrary, the difference 
substantially grows for the remaining parameters. As degeneracy occurs 
using analytical and numerical RFSM5 methods separately, it cannot safely 
be considered as an artifact of the corresponding algorithm. It might be 
either a real feature of the problem of minimization, or a consequence of an 
incomplete domain, towards central regions, of SDH density profiles, owing 
to numerical artifacts (mainly two-body relaxation). The latter possibility 
could be tested by the next generation of high-resolution simulations. 

The current attempt is limited to GPL density profiles, defined by Eq. (1), 
but different alternatives have been exploited in the literature, such as the 
family of density profiles: 

; x = ^ ; (46) 

a 

where the choice (a,/?, 7) = (2,3, 1) corresponds to the Burkert (1995) den- 
sity profile, which resembles the NFW density profile for r ~ 0.02r^jr (e.g. 
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Figure 2: Histogram of the decimal logarithm of the scaled radius, log^mr, 
with regard to numerical (CMV; top) and best fitting (current paper, Tab. 8; 
bottom) RSFM5 methods. Label captions: I-IMP; E-EMP; B-BMP; blank - 
numerical (CMV) RSFM5 methods i.e. same cases as in the upper istogram. 
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Bullock et al. 2001). 

Another possibility is a profile that curves smoothly over to a constant 
density at very small radii (Navarro et al. 2004): 



where the parameter, A, prescribes how fast the density profile turns away 
from a power- law near the centre. The best fit reads (19 sample objects, 
Navarro et al. 2004): A ^ 3ct^ ^ = 0.17216 T 0.021897. 

On the other hand, the RFSM5 method may be extended to any kind of 
density profile. 

5 Conclusion 

Analytical and geometrical properties of generalized power-law (GPL) den- 
sity profiles have been reviewed, and special effort has been devoted to the 
special cases where GPL density profiles reduce to (i) a double power-law 
(DPL), and (ii) a single power- law (SPL). Then GPL density profiles have 
been compared with simulated dark haloes (SDH) density profiles, and non- 
linear least-squares fits have been prescribed, involving five parameters (a 
scaling radius, ro, a scaling density, po, and three exponents, a, (3, 7), which 
specify the fitting profile (RFSM5 method). 

More specifically, the vahdity of a necessary condition for the occurrence 
of an extremum point, has been related to the existence of an intersection 
between three surfaces, 7^ = fk{rQ,a), k = 1,2,3, in a 3-dimension space, 
(0 ro a 7). Using the algorithm has made also establish that the extremum 
point is a fiducial minimum, while the explicit calculation of the Hessian 
determinant has been avoided to gain in simplicity. In absence of a rigorous 
proof, the fiducial minimum could be considered as nothing but a fiducial 
absolute minimum. 

An application has been made to a sample of 17 SDHs on the scale of 
cluster of galaxies, within a flat ACDM cosmological model (Rasia et al. 
2004). In dealing with the averaged SDH density profile (ADP), a virial ra- 
dius, Tyir, equal to the mean over the whole sample, has been assigned, which 
has allowed the calculation of the remaining parameters. The following re- 
sults have been found, (i) A necessary condition for the occurrence of an 




(47) 



38 



cxtrcmum point is satisfied for 8 sample haloes, and is not for the remaining 
9 together with ADP. In the former alternative, an extremum point of mini- 
mum (EMP) may safely exist. In the latter alternative, the occurrence of an 
EMP cannot be excluded, but only a local minimum can be determined, (ii) 
The occurrence of an EMP implies a sum of square residuals which is sys- 
tematically lower than its counterpart deduced by use of numerical RFSM5 
methods. Accordingly, EMPs may safely be thought of as absolute minima. 
With regard to sample haloes where no EMP is detected, the above result 
maintains in 4 cases (including ADP), while the contrary holds for the re- 
maining 5 cases, (iii) The best fit (with no EMP detected) is provided by 
DPL density profiles for 3 sample haloes. In addition, DPL density profiles 
make a rough, but viable approximation in fitting SDH density profiles. The 
contrary holds for SPL density profiles. 

No evident correlation has been found between SDH dynamical state 
(relaxed or merging) and asymptotic inner slope of the logarithmic density 
profile or (for SDH comparable virial masses) scaled radius. 

Mean values and standard deviations of some parameters have been cal- 
culated and, in particular, the decimal logarithm of the scaled radius, ^mr, 
has read < log^^^ >~ 0.59 and (Tsiog^^.^ ~ 0.59, the standard deviation 
exceeding by a factor 3.3-3.4 its counterpart evaluated in an earlier attempt 
using NEW density profiles (Bullock et al. 2001). If a large dispersion still 
maintains for richer samples, in dealing with analytical RSFM5 methods, a 
low dispersion found in A'^-body simulations has seemed to be an artefact, 
due to the assumption of NEW (or any equivalent choice) density profile. 
It has provided additional support to the idea, that NEW density profiles 
may be considered as a convenient way to parametrize SDH density profiles, 
without implying that it necessarily produces the best possible fit (Bullock 
et al. 2001). 

With regard to RESM5 methods formulated in the current paper, the 

exponents of both the best fitting GPL density profile to ADP, (q;,/3, 7) f» 
(0.3,4.5,1.5), and related averages calculated over the whole halo sample, 
(a, (3, 7) ~ (3.0, 3.9, 1.3), have been found to be far from their NEW counter- 
parts, (a, 7) = (1, 3, 1). The last result, together with the large value of the 
standard deviation, agiog^^i^, been interpreted as due to a certain degree 
of degeneracy in fitting GPL to SDH density profiles. If it is a real feature of 
the problem, or it could be reduced by the next generation of high-resolution 
simulations, still has remained an open question. Values of asymptotic inner 
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slope of fitting logarithmic density profiles, have been found to be consistent 
with results from recent high-resolution simulations (Diemand et al. 2004; 
Reed et al. 2004). 
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A On the occurrence of extremum points of 
minimum to the sum of square residuals 

A necessary condition for the occurrence of extremum points of minimum to 
the sum of square residuals, defined by Eq. (34), is expressed by Eqs. (35), or 
more explicitly as: 





(48a) 


U dye 


(48b) 




(48c) 




(48d) 


J2R,—1 = Q ; 




(48e) 



where, for sake of simplicity, it is intended that the residuals and their first 
derivatives are calculated at a fiducial extremum point of minimum, P = 

By use of Eq. (31), the first derivatives of the i-th residual read: 

^ = -7-W-7)U ; (49a) 



42 



dRi 



c 



-1 ; (49b) 



^ = p^ gi + - — -{xi - xc)ti ; (49c) 
op a 

dRi , . 1 / An \ 

^— = (xi - xc) Si ; (49e) 

a 

where both Si and ti depend on Wi = exp]^Q[a(a;j— x^)], according to Eqs. (38a) 
and (38b), respectively. In addition, the vahdity of the relations: 

dRi dRi 

oa a oxc a op a 

may easily be checked. 

The combination of Eqs. (48b) and (49b) yields: 



Y^Ri^Q ; (52) 

i=l 

and the combination of Eqs. (48d), (48e), (50), and (52), yields: 

n 

Y^RiXi = Q ; (53) 

i=l 

on the other hand, multiplying the left and the right-hand side of Eq. (51) 
by Ri and summing over all the residuals, produces: 

~{ ' da ^ a ^dxc 

+ ^^Y.R^-M = —T.R^i^^-^c) (54) 

which, owing to Eqs. (48a), (48c), (48d), (52), and (53), yields: 

Y^x.R^^^O ; (55) 
7^1 9xc 
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similarly, the combination of Eqs. (48a), (49a), (52), and (55), yields: 

n 

Y.R^U = Q ■ (56) 

i=l 

and the combination of Eqs. (48d) and (49d) yields: 

n 

J^RiSi^o ; (57) 
1=1 

finally, the combination of Eqs. (48c), (49c), (56), and (57), yields: 

n 

Y^RiXiti^O ; (58) 
1=1 

in conclusion, Eqs. (52), (53), (56), (57), and (58), make a system of five 
independent equations which is equivalent to (48). 

The substitution of Eqs. (31) and (38a) into (52) yields: 

{y - yc) + 7{-x - xc) + ^ ; (59) 

a 

where overbars denote mean values, according to Eq. (37b). The compari- 
son between alternative expressions of yc, deduced from Eqs. (31) and (59) 
produces: 

Ri^ {yi-y)+l{xi-x)^ {si-s) ; (60) 

a 

where no dependence on yc occurs. 

The substitution of Eq. (60) into (53) yields: 

/3 - 7 

Vs{x, y) + ^Vs{x, x) H Vs{x, s) = ; (61) 

a 

where the function, Vs, is defined by Eqs. (37), and can be conceived as a 
two-dimension (2D) empirical variance. This is why, in the limit of a one- 
dimension (ID) empirical variance, Vs{a^a) = (a^) — (a)^, according to the 
standard notation. It is worth of note that, while ID empirical variances are 
always non negative, the contrary holds for 2D empirical variances. 
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The comparison between alternative expressions of x = {(3 — 'y)/a, de- 
duced from Eqs. (60) and (61), produces: 

Ri = {vi -y)+ - x) 

where no dependence on (3 occurs. Then the combination of Eqs. (56) and 
(62) yields: 

Vs{x, s)vs{y, t) - Vs{s, t)vs{x, y) = 

'j[vs{s, t)vsix, x) - v,{x, s)vs{x, t)] ; (63) 

and the combination of Eqs. (62) and (63) allows an expression of which 
depends only on xc and a, via Eqs. (38a) and (38b). 
The combination of Eqs. (57) and (62) yields: 

Vs{x, s)vs{y, s) - Vs{s, s)vs{x, y) = 

-f[vs{s,s)vs{x,x) - Vs{x,s)vs{x,s)] ; (64) 

which has the same formal expression as Eq. (63), provided the following 
correspondences hold: {y,t) <-> {y,s); {s,t) (s,s); {x,t) <-> {x,s). 
The combination of Eqs. (58) and (62) yields: 

Vs{x, s)vs{y, xt) - Vs{s, xt)vs{x, y) = 

'^[vs{x,x)vs{s,xt) - Vs{x,s)vs{x,xt)] ] (65) 

which has the same formal expression as Eqs. (63), and (64), provided the 
following correspondences hold: {y,t) {y,s) {y,xt); {s,t) (s, s) 
{s,xt); {x,t) ^ {x,s) {x,xt). 

Then we are left with three alternative expressions of 7, which may be 
deduced from Eqs. (63), (64), (65), and are expressed by Eqs. (36a), (36b), 
(36c), respectively. 

A necessary condition for the occurrence of extremum points of minimum 
to the sum of square residuals, defined by Eq. (34), has been formulated via 
Eqs. (48). The condition is also sufficient, provided the Hessian determinant 
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satisfies Eq. (33). The derivation of botli sides of Eq. (35) yields: 



\dujUk)p, 



2E 

1=1 

= 2E 



i=l 



duj J p. 



+ Ri 



duk^ 



+ Ri 



dujduk 



(66a) 
(66b) 



where Ui = xc, U2 = yc, = a, 1*4 = /3, 1*5 = 7. 

The second derivatives of the i-th residuals may be determined, using 
Eqs. (38a), (38b), (38c), (39a), (39b), and (49). The result is: 



d^Rj 
dye 

d^Rj 

d^Rj 
2 

^^R^ 
dycdxc 
d^Rr 



a(/3-7)— InlO ; 

Wi 








; 

/3-7 



(Xi 



a 



X 



^2 

xc)ti + {xi- xcf^lnlO 



dycdl 

d^Rj 
dycdf3 

^^R^ 
dycda 
d^Rj _ 
d'jda 

d'ydxc 



= 
= ; 



= 

= 
1 



Xi - Xc 



a 



tr 



ti_ 



(67a) 
(67b) 
(67c) 
(67d) 

(67e) 
(67f) 
(67g) 
(67h) 
(67i) 
(67j) 
(67k) 
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d(3dxc 



-ti 



1 



Xi - Xc 



ti 



(67m) 



(671) 



df3da 
d'Ri 
dxcda 



{P-^)ix,-xc)-^lnlO ; 



a 



(67n) 



which allows the calculation of the Hessian determinant. The related, ana- 
lytical expression appears to be extremely complicated, and for this reason 
it shall not be calculated here. 

In principle, the above procedure could be used for minimizing the sum 
of absolute value residuals, but at a higher price. This is why the first 
derivatives of the function: 



are not defined at i?j = 0, 1 < i < n, which makes a necessary condition 
for the existence of an absolute minimum, expressed by Eq. (32). satisfied at 
any point, P*{xQ,yQ,a*, where Ri — 0. To avoid this inconvenient, 

our attention has been restricted to the minimization of the sum of square 
residuals. 



n 



F{xc,yc,a,P,^) = J2 \Ri 



(68) 



i=l 
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